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Abstract 

Motivated by the recent experimental observations, we discuss the 
freeze-out properties of the fireball created in central heavy ion colli- 
sions. We find that the freeze-out conditions, like temperature, veloc- 
ity gradient near center of the fireball, are similar for different colliding 
systems and beam energies. This means that the transverse flow is 
stronger in the collisions of heavy nuclei than that of the light ones. 

The system that created in relativistic heavy-ion collision can have both 
longitudinal and transverse expansion (see, e.g., |]]). In order to study the 
hadronic experimental data 0, |3|, one needs a mathematical description of 
a final stage of collisions. At the region near mid-rapidity the boundary 
effect, due to the finite length of the hydrodynamic tube in longitudinal 
direction, can be neglected |4j]. It is possible to use Bjorken's model for 
longitudinal expansion: v\ = xi/t. The same quasi-inert ial flow is inherent 
to Landau model at freeze-out stage || |6|. In this approach, the parameter 
r was introduced to describe proper time of the expanding system. For 
the transverse expansion, we will use a rather general picture proposed in 
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Ref. where, due to the cylindrical symmetry, one has vt{0) = with the 
derivative of the transverse velocity near the center of the fireball v' T (0) ^ 0. 
The transverse velocity increases monotonously as a function of radius r. 

The single particle spectra in a pure hydrodynamic picture without res- 
onance decays are expressed by the integral of the Wigner function over the 
freeze-out surface: 

rl 3 N r 

P°^= da^f(x,p) (1) 



d 3 p 

The Wigner function f(x,p) = fth( x ,p) • p{ x ), where fth{ x ,p) is the lo- 
cal thermal distribution function with a temperature parameter T = 1/13, 
chemical potential fi and 4- velocity field u(x): 

fth(x p) = & + V I ( 2 ) 

The function 



p(x) oc exp 



a 2 
7T ( M ( r > Vl) ~ u(0,y L )) 



exp [—a (cosh?/r(r) 



(3) 



is introduced to describe the finiteness of the relativistically expanding sys- 
tem in transverse direction |7[] . Note that the form of the transverse rapidity 
Dt dependent density is similar to one used in Ref. ||. Here yi and yr are 
the longitudinal and transverse rapidities respectively and a is the intensity 
of transverse flows 



a = K(0) • R T )- 2 . (4) 

For very intensive relativistic transverse flows, a — > 0. On the other hand 
for a>lwe have non-relativistic transverse flows and 

p[x) oc exp(-r /2R 2 T ) (5) 

for yr(r) <C 1. 

Using the saddle point approximation , we have obtained from Eq. (|l|), 
for nixA 3> 1, mr = Jm 2 + p 2 ^ fTj: 
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rf2jV K e -(/3m T+a )(l-4)V2 ^ 



where transverse velocity at the saddle-point r(p) is 

v T = tanh y T (r(p)) = a ^ PT . — (7) 

3m T + a 

Within the model described in Eqs.(l)-(3), the transverse momentum 
distribution (6) is not very sensitive to the details of the velocity profile, see 
Ref. [5]. In the region where 

(3m + a) 2 

m r ~ m < ^— — - (8) 



one has the relationship: 



{3m T + a)(l - 4) 1/2 = (Pm + a) + ~ "° (9) 

(3m + a) 



Therefore we have the simplified Eq. (^j) as: 

(i 2 iV 



mTdmrdy 

witbQ 



e ff(m T -m) 



(10) 



T7? 

l/& ff = T eff = T+-. (11) 

a 

The latter equation connects the slope of the transverse mass spectra with 
the intensity of transverse flows 

jy 2 hydrodynamical length 
transverse radius 

The intensity of transverse flows a is defined by the Eq. (|4]) and does not 
depend on m, it depends only on A, where A is the nucleus atomic number. 



1 The linear dependence of the slope parameter of the transverse momentum distri- 
butions on particle mass was obtained early for spherically symmetric non-relativistic 
expanding system in Ref. || . 
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The system finally decays when the rarefaction wave reaches the center 
part of the fireball. The analysis of the experimentally measured slope pa- 
rameters (T e ff) as a function of particle mass ]TU| with two parameters, T 



and a, shows that the freeze-out temperature T is approximately constant 
for different colliding systems at the beam energy > 10A-GeV. Hence it is 
natural to conclude that the physical condition near the center is the same 
for different colliding systems, i.e. v' T (0) is constant. From Eq.([|) we obtain 

V ' T (0) = (cdf T y 1/2 (13) 

In our approach Rt is the Gaussian-like radius of a decaying system. We 
suppose that Rt depends on A only. It means that in the same collisions 
a freeze-out "size" is unique for different particle species. The value Rt 
is connected with the so called sideward radius R s |TT|, obtained from the 



fit of the experimentally measured two-particle correlation functions, in the 
following way. When v T = pj^ +a "C 1 (non-relativistic approximation) it 
can be shown that only v' T {0) has an influence on the interfere metry radii 
and any details of the transverse velocity profile Vt(t) = tanh^T(^) are 
not important |7|, |12|| . Within this approximation, one has the following 
expression for the Gaussian transversal sideward radius |7], [12| : 



From this equation we evaluate Rt, using experimental data for Rs(pt), 
and (3 and a from the fitting procedure [|l(J based on Eq. (|TT]) . To min- 
imize the influence of the resonance decays, the interferometry radius for 
every analyzed particle species has to be measured at the point where pt 
is large enough. On the other hand, in order to provide a validity of the 
non-relativistic approximation and the correctness of the condition Eq. (|8|), 
the value of pt should be limited with v T <C 1. 

Using the experimental data, listed in Table I, we evaluate the freeze-out 
temperature and velocity gradient at the center of the fireball. The results 
of the fit are shown in Fig. 1. Note that there are two important features 

2 Assuming the validity of this approach, we found that, even for the heaviest colliding 
system Pb+Pb [10], a = 6.57 and vt — 0.32. These values were obtained at tut ~ 0.5 
GeV where is the measured highest pair momentum region. For smaller pt and lighter 
colliding system, this assumption works better. 
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in the plot. First, the values of the intrinsic temperature T seems to be a 
constant (top plot of Fig. 1). It does not dependent on the size of colliding 
system. Furthermore, for collisions at the AGS energies (11 - 15 A-GeV/c) 
and SPS energies (158 - 200 A-GeV/c), this parameter is approximately the 
same. Secondly, the value v T (0) is also closed to a constant (see bottom plot 
of Fig. 1). 

It has been noted that the temperature is limited for collisions at beam 
energies larger than 11 A-GeV (Fig. 2 of Ref. |H|). One should note that, 
although we do not observe the difference in v T (0) between the collisions at 
AGS and SPS energies, it is quite possible that at lower beam energies one 
gets somewhat different values. The saturation of the freeze-out temperature 
might reflect the limiting temperature hypothesis pointed by Hagedorn some 



years ago [15], |16| . On the other hand, the constant behavior of the velocity 
gradient near the center of the fireball is a new result yield by the present 
study. It clearly shows that at the final stage, the freeze-out conditions are 
approximately the same for different size and different beam energies of the 
collisions. Constant transverse velocity gradient actually means that both 
the averaged transverse velocity and flow intensity grow with the size of the 
colliding nuclei. Indeed, the average transverse hydrodynamic velocity that 
can be calculated using the final distribution function (H) is 

= J^Jda^v T (x)f(x,p) = Jda,u^x)p(x)v T (x) 
f^fda^f(x,p) " Jda^(x)p(x) 

For r(r) = const at the freeze-out hypersurface we can get from Eq. (|15|) : 

^ ftanhy T (r) coshy T (r) exp [-a cosh y T {r)} rdr 
J coshyT(r) exp [—a cosh|/r(r)] rdr 

Using saddle point approximation, in the limit of a ^> 1, one can show that 
the average transverse hydrodynamic velocity is 

. . , . . f exr>(—r 2 /2Kl)r 2 dr 



/ exp(— r 2 /2R T )rdr 



(17) 
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Let us note that vt 7^ {vt)- The physical meaning is that, in the saddle 
point approximation, vt is the transverse velocity of the fluid element which 
gives the main contribution to transverse momentum spectrum at given pj*. 
On the other hand, [vj) is the averaged transverse hydrodynamic velocity 
that characterizes whole system. From Eqs. (16) and (17), one can see 
that a is directly connected with (i>t)- So we find that, in approximation of 
Eq. (17), the averaged values of the transverse velocity (vt) are 0.49c and 
0.33c for Pb+Pb and S+S central collisions [10], respectively. Formally the 
averaged velocity (vt) is consistent with zero for p+p collisions although the 
application of the thermal model for such collisions remains an open question. 
At the freeze-out surface, the effective slopes of spectra for different colliding 
systems and different mass particle species can be presented by the non- 
relativistic average transverse hydrodynamic velocity: 

T eff ^T+-m(v T ) 2 . (18) 

7T 

The merit of this equation is that by fitting the measured slope parameter 
as a function of particle mass, one can separate the collective motion from 
the thermal motion. Therefore, the intrinsic freeze-out temperature T and 
averaged collective velocity (vt) can be readily extracted from the data. The 
measured slope parameters of pions, kaons, and protons [10] seem to obey 
the linear mass dependence as given by Eq. (18). One new finding from 
this study is the linear tie of the averaged transverse velocity (vt) with the 
radius Rt, see Eq. (17). It is worth mentioning that the hydrodynamic 
arguments for the behavior of the slope parameter as a function of a particle 
mass is valid as long as the considered particle species remain a part of the 
fireball and hence are participating in the frequent rescatterings. For those 
particles with high probability of destruction, such description is not valid 
anymore. It is possible that those particles freeze-out earlier than the other 
hadrons which participate in the evolution of the system longer. They do 
not have enough time to acquire a common collective velocity and hence the 
corresponding slope parameters could be smaller than the predicted by Eq. 
(18). Indeed, the recent preliminary results reported by WA97 |17| and NA49 



[ |I~8|j show such deviations for hyperons, particularly for Q particle ||19[| . In 
order to understand the dynamics involved for those particles, one has to 
take the collision rate for individual considered particle into account. 

In summary, using a thermal model we analyzed the recent experimental 
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data from heavy ion collisions. We found that the freeze-out temperature 
is a constant for all colliding systems at Eb eam >10 A-GeV and a constant 
velocity gradient at the center of the fireball (r — > 0). 
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System 


1/v 7 ^ 


Rs 


P + P [10] 


0.01±0.11 


0.8±0.2 


S + S [10] 


0.26±0.10 


2.9±0.2 


Pb + Pb [10] 
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Si + Al [13] 


0.29±0.12 


2.56±0.17 


Si + Au [13] 


0.30±0.12 


2.4±1.1 


Au + Au [14] 


0.42±0.11 


3.57±0.52 
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Figure 1: The freeze-out temperature T (top) and flow velocity gradient 
(bottom) calculated at r « 0. The set of data for this calculation is the 
taken from the references listed in Table I. 
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